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We study the process of low-energy electron capture by the SF6 molecule. Our approach is based 
on the model of Gauyacq and Herzenberg [J. Phys. B 17, 1155 (1984)] in which the electron motion 
is coupled to the fully symmetric vibrational mode through a weakly bound or virtual s state. By 
tuning the two free parameters of the model, we achieve an accurate description of the measured 
electron attachment cross section and good agreement with vibrational excitation cross sections of 
the fully symmetric mode. An extension of the model provides a limit on the characteristic time 
of intramolecular vibrational relaxation in highly-excited SF^" . By evaluating the total vibrational 
spectrum density of SF^~, we estimate the widths of the vibrational Feshbach resonances of the 
long-lived negative ion. We also analyse the possible distribution of the widths and its effect on the 
lifetime measurements, and investigate nonexponential decay features in metastable SF^". 
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I. INTRODUCTION 



Electron attachment to SFg is a fascinating problem. In spite of a lot of attention and a wealth of experimental 
I | data some basic questions, e.g. that of lifetimes of metastable SFg , lack definitive answers. In this paper we show 
that the electron capture process is described well by a zero-range-type model [2]. We determine the parameters of 
the model by comparing the theory with experimental data on the attachment, total and vibrational excitation cross 
sections. We then study the autodetachment widths of SFg~, and analyze its lifetimes and nonexponential decay. Here 
the experimental situation is less clear. Our calculation yields characteristic lifetimes of about a millisecond, using 
possibly the most accurate set of SFg~ vibrational frequencies Q . We investigate the nature of nonexponential decay 
• i-h ' of metastable anions due to level-to-level fluctuations of the widths and a distribution of the incident electron and 
target energies. 

Sulfur hexafluoride (SFg) has long been know as an electron scavenger because of its large low-energy electron 
p> ' attachment cross section. This feature of SFg is important for its applications as a gaseous dielectric and makes for 
the rich physics of low-energy electron collisions with it. The energy dependence of the electron capture cross section 
is well established experimentally P, d, H, H, 0, B @- Below 10 meV it shows 1/v behavior characteristic of s-wave 
inelastic processes. At higher energies towards 100 meV it approaches 60% of the unitary limit ir/k 2 for the reaction 
' cross sections, where k is the incident electron momentum (atomic units are used throughout). In addition, SFg also 
has a large elastic scattering cross section reaching ~ 10 3 a.u. at electron energies of a few meV, which can be inferred 
7— I ■ from the measured total cross secti on IToL [Tl| . Note though that experimental data on low-energy elastic collisions, 
including differential cross sections [12|, and inelastic scattering cross sections, e.g. those of vibrational excitations 
£NJ ■ [1, EH E3] , are relatively scarce. 

Low-energy electron attachment leads primarily to the formation of long-lived parent negative ions, 

oo : 

O e-+SFg^SFg-*, (1) 



At electron energies e > 0.2 eV (and below 3 eV) dissociative attachment, 



e- + SFg -> SF 5 - + F, (2) 
C3 ' 

becomes dominant. Although this process is usually regarded as a channel separate from JT]), some recent evidence 
suggests that SF^~ can be formed in the decay of metastable SFg [HI, El] , 

SFg * -> SFg + F. (3) 
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However, at low energies the lifetimes of SF 6 * in the absence of collisions are limited by electron autodetachment, 



SF 6 -* ^SF 6 + e-. (4) 

Numerous SFg"* lifetime measurements fPA [l8l [l9l 21^ 22 ^ 28l 24 ,_ 25l 26j show considerable variation, depending 
on the technique used. Time-of- flight experiments 17 TTTg . 19L l20l 2lL |22(| typically yield values of several tens of 



microseconds, while ion-cyclotron-resonance methods give values of about a millisecond or larger 

Using the ion-cyclotron-resonance method, Odom et al. [24] found that the apparent lifetime of SFg"* varied as a 
function of the observation time, and surmised that SFJ were formed in a distribution of states with different lifetimes. 
The most recent study that used a Penning ion trap [26| also points to the formation of ions with a range of 1-10 ms 
lifetimes. A similar conclusion was derived from the time-of-flight measurements by Delmore and Appelhans [13, Hl| 
at microsecond time scales. Their analysis of the SFg"* decay indicated multiple lifetimes or groups of lifetimes in 
the interval between 2 to 30 /is, and showed that the population of states with different lifetimes depended on the 
temperature of SFg molecules. On the other hand, Ref. [22| where a free jet expansion was used to cool down the SFg 
molecules, reported a single lifetime of 19.1 ± 2.7 /is. All these differences are usually attributed to the differences in 
the experimental conditions under which SFg" * are formed, i.e., the incident electron energy and the internal energy 
of the target molecule [l|. However, detailed understanding is still lacking. 

In contrast to the problem of lifetimes, the process of low-energy electron attachment to SFg is described well 
by a simple zero-range- type model of Gauyacq and Herzenberg jjj. According to this model, the incoming s-wave 
electron undergoes strong resonant scattering on a virtual or weakly bound level of the SFg molecule, with a near- 
zero energy. Electron trapping occurs via population of this fully symmetric state which is strongly coupled to the 
symmetric stretch ("breathing") vibrational mode v\. Electron capture initiates the motion of the fluorine nuclei 
towards the equilibrium configuration of the negative ion. This process is accompanied by rapid intramolecular 
vibrational redistribution (IVR) of the breathing mode energy among other vibrational degrees of freedom. As a 
result, the probability for the nuclei to return to the equilibrium configuration of the neutral SFg becomes small, and 
long-lived metastable anions are formed. 

In this paper we perform a comprehensive study of electron attachment using the zero-range model. The attachment 
cross section is sensitive to the behavior of the SFg and SFg" potential energy curves, as a function of the S-F bond 
length, near their merging point. Details of this behavior are incorporated in the model via two parameters, namely 
the energy of the virtual (or weakly bound) level and its coupling to the breathing mode. We develop a new effective 
matching procedure that connects the region of electron capture near the merging point, where the zero-range model 
can be applied, with the outer region of adiabatic semiclassical nuclear motion. We also generalize the model to 
include the possibility of the nuclear framework to oscillate back to its initial configuration, in order to study the 
influence of the rate of IVR on + SFg collisions. Our aim here is to test the model by comparison with experimental 
data on the attachment, total and vibrational excitation cross sections, and thus determine its parameters. Another 
goal is to compare these parameters with the results of potential curve calculations. The SFg and SFT potential 
curves have been established quite well by now [HH?]], overcoming earlier uncertainties [28L l29l l30l. l3~il. l32f . 

A somewhat different theoretical approach was taken recently by Fabrikant et al. [9]. It starts from the R- 
matrix formalism, and goes beyond the model of Gauyacq and Herzenberg by including the long-range polarization 
potential — a/2r 4 and dipole coupling between the s and p waves. Contributions of higher electron partial waves 
are also included, to describe the elastic and vibrational excitation cross sections more accurately. In this theory 
the two parameters of the s-wave coupling to the V\ mode are allowed to be complex. Their values are found by 
fitting the experimental attachment and total cross sections. However, the use of complex parameters makes them 
phenomenological, and the physical connection with the molecular potential curves is lost. 

The second part of our paper concerns the evaluation of lifetimes of SFg~* due to electron autodetachment. The first 
estimates of the metastable anion lifetimes and their relation to the attachment cross section, vibrational spectrum 
density of SFg" and electron affinity of SF 6 were made in Ref. [l8j]. The rate constants of the processes © and 
(HI) , were later studied [H, [34], [35|, |36| using the quasi-equilibrium (or RRKM) theory (3?j]- The approach of Refs. 
fl8L l33l . l34| is based on the principle of detailed balance, and requires the knowledge of the attachment cross section 
and the anion vibrational spectrum density. The standard RRKM approach is similar, assuming in addition that the 
transition probability is unity, i.e., that the attachment cross section is equal to its unitary limit. On the other hand, 
the RRKM requires the knowledge of the density of so-called transition states. The autodetachment lifetime depends 
strongly on such parameters as the electron affinity E a and vibrational spectrum density of SFg", which were not 
known well. Using different sets of data, lifetimes from microseconds to milliseconds were obtained in a wide range 
of incident electron energies. It should be noted that these methods yield the detachment rates (inverse lifetimes) 
averaged over the ensemble of metastable anion states. The distribution of lifetimes within such ensemble that may 
cause a variation in the observed lifetimes, has not been studied. 

In the present paper we analyze the dependence of the autodetachment rate on the incident electron energy in 
the interval from zero to 300 meV for different target temperatures. Our calculations are based on the accurate 



values of the attachment cross section, and take into account the distribution of the target SFq molecules over the 
rotational and vibrational states, as well as the rotational and vibrational degrees of freedom of the SF^ anion. We 
also analyse fluctuations of the decay rate over the ensemble of metastable SF^ , to see if they can explain observations 
of nonexponential decay. 



II. ELECTRON SCATTERING AND CAPTURE 



A. Attachment model 

Following the approach of Gauyacq and Herzenberg Q we employ the zero-range potential (ZRP) model to describe 
the electron interaction with SFk. ZRP theory is a well-known tool suitable for problems of low-energy electron 
scattering and negative ions [33, [39| , especially in those cases when the cross section is enhanced by the existence of 
a shallow bound state or low- lying virtual level with zero angular momentum (i.e., in the s wave). Application of the 
ZRP method to the e~ + SF 6 system provides the wavefunction for the continuous spectrum electron and SF 6 in the 
vicinity of its equilibrium configuration. On the other hand, the wavefunction of SFg~ formed as a result of electron 
attachment can be described in the Born-Oppcnhcimcr approximation, treating the nuclear motion semiclassically. 
Matching of these wavefunctions yields the electron attachment cross section together with those of elastic scattering 
and vibrational excitation. 

The electronic state responsible for the capture process is a fully symmetric weakly bound state of + SF6. Outside 
the molecule the electron wavefunction ipo is spherically symmetric, and is given by 



V Z7T r 

where r is the electron coordinate and k is related to the bound state energy £o = — k 2 /2. The energy £q depends on 
all nuclear coordinates. However, due to its symmetry, this electron state is most strongly coupled to the breathing 
(symmetric stretch) vibrational mode of SFg. This allows one to consider Sq and k as functions of the normal coordinate 
q of the breathing mode, q = R — i?o, where R is the S-F bond length, and Rq is its value at the equilibrium of SFg. 
Near the equilibrium, n(q) can be expanded in a power series. Keeping the first two terms [2j, 

K(q) = K Q + K 1 q, (6) 

where m > 0, since the binding increases with the increase in R (R = 3.25 a.u. at the equilibrium of SF^", while 
Ro = 2.96 a.u. Q). 

The negative ion bound state ([5|) exists when n > 0, i.e., for q > qo = —kq/ki, while k < corresponds to a virtual 
state [40|. The absolute value of k is expected to be small, |«o| "C 1 a.u., because for the electron capture to be 
effective, the scattering length k^ 1 should be large in the Frank-Condon region. The sign of kq indicates whether the 
negative ion state is real or virtual at the equilibrium of the neutral. The two cases are illustrated in Fig. [TJ which 
shows the potential curves of SFg and SF^~ for the symmetric stretch coordinate. 

To the right of the merging point, qo, indicated in Fig. Q] by vertical arrows, the negative ion energy Uo(q) is 
given by the sum of that of the neutral SF 6 and the bound state energy £o(q) = —n 2 {q)/2. The SF 6 potential curve 
shown was calculated in Ref. [13]. Near the equilibrium it can be approximated by ^Mw 2 g 2 , where uj and M are 
the frequency and mass of the breathing mode. Using Eq. §E§ one obtains a quadratic approximation for the SFg~ 
potential curve near the merging point, 

Uo{q) . M^! _ 0^0 + ^, (7) 

For q to the far right of the merging point, Uo(q) should approach the anion potential energy curve, e.g., that calculated 
in Ref. 

In this model the electron scattering length k -1 depends on the breathing vibrational coordinate q. This means that 
the electron interacts only with this particular mode, and no other vibrations are excited in the process of electron 
scattering or capture. Hence, all nuclear coordinates except q can be omitted, and the the total wavefunction of the 
system can be written as 

*(r, g) = e^ Xna (q) + £ ye ifc ^ X „(g). (8) 




Coordinate q (a.u.) 

FIG. 1: SF6 and SF^" potential curves for the symmetric stretch coordinate. Dashed and dotted curves are the SF6 and SFg~ 
energies, respectively, calculated using the 2nd order M0ller-Plesset perturbation theory [27|. Dot-dashed curves show the SF^ 
energy in the quadratic approximation for go = —0.04, fti = 2.0 (a) and qo = 0.034, ki =4.1 (b). Solid curves interpolate 
between the quadratic and numerical SF^" potential curves and represent the lowest adiabatic energy of the e~ + SFg system. 
Horizontal dashed lines show the lowest total energy of e~ + SFq collision, uj/2. 

The wavefunction ([5]) is an expansion over the SFg breathing mode vibrational states Xn(l), taken in the harmonic 
approximation, with energies E n = o;(n+g). The first term in Eq.([H]) describes the electron with momentum k = k no 
incident on the target in the initial state no- The sum over n accounts for elastic scattering (n — no) and vibrationally 
inelastic processes (n ^ no). Energy conservation, ifc 2 + E no — |fc 2 + E n , determines the corresponding electron 

momenta, k n = \/k 2 — 1uj(n — no). 

Note that the sum in Eq. ([5]) includes both open (real k n ) and closed (k n = i\k n \) channels. In the former, 
the electron escapes but the nuclear motion is finite. Closed channels involve the electron localised near the origin, 
with the vibrational motion in progressively higher n states. Its contribution describes electron attachment, with 
the nuclei sliding down the negative ion potential curve towards greater q. It takes the form of the SF^ anion 
adiabatic wavefunction \E , at t(r,g) (see below). Of course, the true anion potential curve (Fig. [T|) does not allow for 
the infinite nuclear motion, as the nuclear framework swings back to the neutral equilibrium after one vibrational 
period. However, in a molecule with many vibrational degrees of freedom the energy may dissipate into other modes, 
providing for electron capture on much longer time scales. 

The electronic parts of the wavefunction are plane or spherical waves, which ensures its correct asymptotic form. 



They are solutions of the free-particle Schrodinger equation. This is in accordance with the ZRP method, in which 
the potential affects the wavefunction through the boundary condition at the origin, 



1 9(rV) 

sV „ - -"<* < 9 » 

Because of the r — > limit, the ZRP affects only the electron s wave. 

Applying Eq. ([9|) with n{q) from ([6]) to the wavefunction {5J and projecting the resulting equation onto each of the 
nuclear vibrational states Xm one obtains a set of linear equations for the amplitudes /„ (n — 0, 1 , . . . ) 



(ik n + K )fn + i , = (Vnfn-1 + VnTTfn+l) = -5 nno . (10) 

The general solution of this 2nd-order recurrence relation is a linear combination of two independent solutions with 
arbitrary coefficients. One "boundary condition" is provided by Eq. (|10|) with n = 0. The other boundary condition 
is set at large n. It is related to the asymptotic behaviour of 4 , a tt(r, q) at large q. Analysis presented in Appendix lAl 
shows that there is a simple relation between the nuclear coordinate q and the quantum number n of the terms in the 
expansion, Eq. l[8|). which contribute significantly to the wavefunction at this q, 

n&n + [k 2 /2-£ (q)}/cj. (11) 

Hence, the wavefunction \&(r, q) from Eq. (JS} must be extended to the region where the incident electron is bound, 
and matched at some point q = g m with the negative ion wavefunction ^ r a tt(r, q). The choice of q m should not affect 
the capture cross section. Physically, it is restricted to the range of validity of expansion (JBJ), so q m should not be too 
large, but sufficient to neglect nonadiabatic effects. 

In this region one can write \P a tt(r, g) in the Born-Oppcnheimer approximation, 

*^(r,5)=^(r,g)x(g), (12) 

where ipo( r ,q) is the bound electron wavefunction (|5"j). which depends on the nuclear coordinate via k — n(q), Eq. 
©, and x(q) is the wavefunction of the nuclear motion of the anion. It is in general a superposition of the outgoing 
and incoming (reflected) waves, 

X(q)=A X M (q)+B X (- ) (q), (13) 

where X can t> e written explicitly in the semiclassical (WKB) approximation 0, as x (±) = u-V2 e ±i/pd« 5 p and 
v being the classical nuclear motion momentum and velocity. 

The amplitude of reflection, R = B/A, depends on the behavior of the negative ion term Uo(q) far from the merging 
point and also on the coupling of the anion breathing mode to the other vibrational modes (i.e., on the effectiveness 
of IVR). In the absence of such coupling one has \R\ = 1 and no capture takes place. In treating electron attachment 
to SFg, the usual assumption is that this coupling is strong, i.e. \R\ -C 1, so that only the outgoing wave is retained 
in Eq. (HU) 0). 

In order to find the solutions of Eq. (jTUJ) corresponding to X ^(q), note that for a sufficiently large n one has: 
ik n + kq ~ —\2um, and Eq. (|10|) turns into a recurrence relation with constant coefficients. Hence, asymptotically, 
the two independent solutions of Eq. (JTUJ) behave as f n oc £™, where £ is either of the two complex-conjugate roots of 
the corresponding characteristic equation. It turns out (see Appendix |X| that one of these solutions corresponds to 
X ^ + \ while the other to 

In practice, Eq. (fTQ|) approaches the constant-coefficient limit slowly, as the coefficients vary slowly with n (4ll |. 
However, this slowness does allow us to truncate the set (fTU|) at some large n — N by using /;v+i = £/jv, where £ 
satisfies the quadratic equation, 



N+l+ 2 n , , s _ . IN 



Ki\ —C -(\kN\-K )C + Ki\ =0, (14) 



and where — iy // 2uj(N — no) — k 2 = i\hjf\ is imaginary. The two roots of Eq. (I14|) are 

(15) 



VMuj(\k N \ - k ) . . 

C± — — = =■= % 

^/2{N+l) Kl 



N Muj(\k N \ - k ) 2 



N + l 2{N+\)k\ 



The choice of N is related to the matching point q m by means of Eq. (fTTj) . The expression in square brackets is 
real because the anion potential jJJ must decrease at large q, which requires n\ > Alui 2 . Note also that |£±| 2 = 
•v/ N/ (N + 1) ~ 1 for large N, which means that £± ~ e ±l13 is just a phase factor. 

As follows from Appendix [A] solving the set of N + 1 equations l|10p (n = 0, 1, ... , N) with the additional 
condition fy+i = £+/iV> generates a set of amplitudes (we denote them fn) that corresponds to the part of 

the nuclear wavefunction (| 1 3|) . On the other hand, using /jv+i = £— /jv> on< 3 obtains a set of amplitudes, /n , which 
corresponds to B%y~\ There is a simple linear relation between the amplitudes and /n \ and i? at large n: 
ft ] /fr[ +) =Rf = e- i5n R, where S n is a phase, see Eq. (|A"20t . 

In the general case, 0, the solution / n of the recurrence relations (jTUJ) that matches \& a tt(r, 5), Eqs. (fT"2"|) . (1 13[) . 
has the form 

/n = C /(")), (16) 

where C is a normalization constant. Therefore, we can set /jv+i = C(l + i?/) and /at = C(^ 1 + Rf^Z 1 )- The 
amplitudes / n for n = N — 1, . . . , no are then found successively from Eq. (jTUJ) . The same is done for f n with 
n = 0, . . . , no, starting from an arbitrary /o, e.g., fa = 1. The two solutions are then matched at n = no, and the 
overall normalization is determined by substitution in the inhomogeneous equation (|10| with n = no- 

In practice, Eq. (jTUJ) can be truncated at relatively low n = N, e.g., AT ~ 10 can be used for electron energies below 
350 meV without significant error. By means of Eq. (|11|) . this allows us to keep a physically meaningful value of the 
matching coordinate, q nl — 0.08, see Fig. Q] 

The elastic (n = no) and vibrationally inelastic (n 7^ no) cross sections are given by 

47r 2 

c„ = — fc n \f n \ , (17) 
where k n is real for the open channels. The total cross section is given by the optical theorem [4(1 , as 

47T 

0t = — Im/„ , (18) 

and the attachment cross section is obtained as 

o-att = y (im /„ - ^ |/„| 2 Re fc„) . (19) 

n 

It can also be found directly from the asymptotic behaviour of /„ for closed channels, Eq. (|A15|) . 

Note that in comparison with Ref. 0, our method allows one to account for nonzero reflection, i.e., for incomplete 
vibrational relaxation, and provides the cross sections of all processes (not just attachment). 



B. Numerical results 



The model described in Sec. Ill Al contains five parameters. The frequency of the SFg breathing mode is well 
established experimentally, w = 96 meV = 3.5 x 10 -3 a.u. [l|, and the corresponding mass is M = 6top = 2.1 x 10 5 a.u. 
The other parameters of the model, namely kq and Ki, which characterise the anion potential curve, and the reflection 
amplitude R, are not known a priori. 

To determine their values and to verify the model itself, we have performed numerical calculations of the cross 
sections in a wide range of parameters, and compared the results with the experimental data on attachment 0,0,11], 
total scattering [ll|, and vibrational excitations @. Our calculations have been done for the target molecules in the 
ground vibrational state of the breathing mode (i.e., no = 0), since at room temperatures the fraction of excited states 
of this mode is small. 

The ZRP model is expected to work best at low electron energies. Here the attachment data is the most accurate 
of all the measured SFg cross sections, and we use it as our main guide in the search for the optimal parameters. 
To characterise the discrepancy between the theory and experiment we calculate the "mean-squared relative error" , 
X] = X)i'=i[^ cr att(£i)/ <7 att(ei)] 2 /36, where Acr at t is the difference between the theoretical and experimental cross 
sections, using 36 energies £{ between 0.1 and 160 meV, as given in the tables of recommended cross sections in Ref. 
[l[. Let us first examine the results obtained for R = (rapid IVR), and then look at the effect of the reflected wave, 
R ^ (incomplete IVR). 




1. Rapid IVR (R = 0) 

The values of r] for R = are shown on a density plot in Fig. as a function of n± and go = —kq/ki. The range 
of parameters is limited by k% > yMu) = 1.6 a.u., and the use of go instead of kq emphasizes the sensitivity of the 
attachment cross section to the position of the merging point of the neutral and anion potential curves. 

Figure [5] shows that the absolute minimum of rj is achieved for go ~ —0.04 and k\ sa 2. A negative value of go means 
that the electron does form a weakly bound SFg~ state at the equilibrium of the neutral, g = 0. Figure [TJ a) shows the 
corresponding anion potential curve, Uo(q) from Eq. ([7]), by the dot-dashed curve. For comparison we also show the 
SFg potential curve calculated in [27| . This potential curve does not reproduce Uq (g) near the merging point [42| but 
should be reasonably accurate at larger q. Figure [T] also shows that the analytical and numerical potential curves can 
be matched. The position of the matching point g m = 0.08 used in our calculations is indicated by an arrow. Note 
though, that for R — the calculation does not require any information about the anion potential curve far from the 
merging point, and the choice of the matching point (or the truncation number N) is not critical. 

The attachment cross section calculated for go = —0.04, ki = 2.0 and R — (dashed curve in Fig. [3]) reproduces 
both the smooth decrease of the measured cross section below the v\ vibrational threshold, e < lu [43j, and its rapid 
drop for e > lo. However, the total cross section calculated with these parameters is noticeably higher than experiment 
for e > u) (see Fig. 2]). 

The other possibility suggested by Fig. [5] is that go > 0, where rj has a second local minimum in the form of a 
narrow "valley". It corresponds to a virtual electron level at the equilibrium of SFq. Here the quality of the fit is not 
very sensitive to the precise value of K\. Choosing go — 0.034, k\ = 4.1 and R = (thick solid curves in Figs. [3] and 
2]) we obtain a much better description of the total cross section, while the fit of the attachment data is only slightly 
worse. In Fig. |Tfb) we show the corresponding behavior of the SFg~ potential curve near the equilibrium. In this case 
the quadratic curve matches the SFg curve from Ref. [13] even closer at q « g m = 0.08. Note though, that due to 
the proximity of g m to go, the adiabatic approximation for the nuclear motion is not as accurate here as in the go < 
case (the required truncation N being lower). This means that the actual shape of the SFg potential curve in the 
vicinity of g m may have a small effect on the cross section. 

Generally, the ZRP model is expected to provide a more accurate description of the attachment cross section than 
the total cross section. Due to the symmetry of the negative ion state and the role played by the symmetric stretch 
mode, the low-energy attachment model can be restricted to the electron s wave. However, contributions of higher 
partial waves to the total cross section (in particular, due to excitation of the strong infrared-active v-$ mode), may 
become sizeable even at low energies. This is indicated by the observed anisotropy of the total cross section [1JJ 
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FIG. 3: Electron attachment cross sections for SFs. Solid squares are the experimental data for metastable SF 6 obtained by 
Hotop et al. as given in Ref. 0]; thin solids curves show the data for SF^T and SF^ obtained at the SFg nozzle temperature of 
300 K Q. Calculations: q = -0.04, m = 2.0, R = (dashed curve); q = 0.034, Ki = 4.1, R = (thick solid curve), R ^ 0, 
7 = 0.1ui a (dotted curve), and 7 = 0.01w a (dot-dashed curve). 
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FIG. 4: Total cross section for electron scattering from SFg. Experimental data: solid squares, Ref. [TT| . Calculations: 
go = —0.04, ki = 2.0, R = (dashed curve); qo = 0.034, ki = 4.1, R = (thick solid curve), R 7^ 0, 7 = 0.1oj a (dotted curve), 
and 7 = 0.01o; a (dot-dashed curve). 

(see also Ref. [9|). Since the present ZRP model does not take into account any of these contributions, one could 
expect that the cross section (fTB|) would be a lower bound for the true total cross section. Hence, we believe that the 
parameters qo = 0.034 and k,\ =4.1 are more realistic, in spite of a less accurate fit of the attachment cross section 
than that for q = -0.04, K\ = 2.0 . 

As a final check of the model for R = 0, in Fig. [5] we compare the calculated cross sections with the v\ and 2v\ 
vibrational excitation differential cross sections measured at the scattering angles of 30° and 135° [9(. Since the ZRP 
model cross sections, Eq. (|T7|) . are isotropic, the differential cross section are found as da n /dil — a n j (47r). As in Figs. 




FIG. 5: Differential vibrational excitation cross sections of the symmetric stretch mode of SFg (top, v\\ bottom, 2l>\). Thin 
solid lines are the experimental results of Ref. || obtained at 30° and 135°. Dashed curves are the cross sections calculated 
for go = —0.04 and ki = 2.0, and thick solid curves are those for go = 0.034 and k\ =4.1 (both for R — 0). 



[3] and HJ the two sets of theory curves in Fig. [5] correspond to go = —0.04, K\ = 2.0, and go = 0.034, K\ = 4.1. The 
experimental cross sections for the two angles are very different close to the threshold. However, the data within 30 
meV of the threshold may not be reliable [4J] , while beyond this region the anisotropy of the differential cross section 
is greatly reduced. Here the experiment is clearly in much better agreement with the calculation for a positive value of 
go (i.e., that for which the anion state is virtual at the equilibrium of the neutral). This is especially clear for the n = 2 
excitation. The calculation for go = 0.034 and «i =4.1 also shows the same cusps at higher vibrational thresholds as 
the experiment. We thus conclude that the experimental data for the low-energy attachment and scattering favour 
the potential curve scheme shown in Fig. [Hb). 



2. Incomplete IVR (R ^ 0) 



In the calculations above we used R = 0, i.e., we neglected the reflected wave ^ m the anion wavefunction (|13jl . 
In general, the size of the reflection amplitude R is determined by the rate of IVR, which depends on the coupling 



between the vibrational modes. A full calculation of the vibrational dynamics of highly excited SFg~ is a nontrivial 
task well beyond the scope of the present work. Hence, we introduce the IVR rate phenomenologically as a width 7 
of the breathing mode. This is done by adding a small imaginary part —ij/2 to the anion potential energy Uo(q). 
The reflection amplitude R is then obtained semiclassically as outlined in Appendix lAl 

In the numerical calculations we treat the anion potential curve for q > q m in the harmonic approximation. The 
IVR damping of the reflected wave is then proportional to ^j/uj a , Eq. (|A19|) . where uj a is the anion breathing mode 
frequency. For strong damping, e.g., j/i-u a — 1, the results are very close those obtained with R = 0. The reflected 
wave here is suppressed by a factor of about w 0.05. In contrast, smaller IVR rates lead to drastic changes in the 
cross sections. In Figs. [3] and H] we show the cross sections for go = 0.034, n\ — 4.1, and two values of the damping 
parameter: -f/oj a = 0.1 and 0.01 [43]. (The effect of damping on the cross sections for qo < is broadly similar.) 

We see that allowing for a sizeable reflected wave results in the emergence of anion vibrational Feshbach resonances 
and overall suppression of the attachment cross section. Both effects spoil the good agreement with experiment 
achieved in the calculations with R — 0. We thus conclude that the IVR in SFg" is very rapid, 7 > uj a . It takes place 
over the time of one vibrational period of the breathing mode, that indicates strong coupling between the breathing 
mode and other nuclear degrees of freedom in SFg . 

In the past, the IVR in highly vibrationally excited SFg molecules was studied in multiphoton laser excitation 
experiments (see, e.g., and references therein). These and other studies [47j indicate that the IVR rate of the 
strong infrared-active 1/3 mode is noticeably smaller than our estimate. However, the two situations are quite different. 
In the capture process the amount of energy equal to the electron affinity of SFg, E a — 1.06 eV (recommended value 
(H), is instantaneously deposited into a single mode. This mode then has a much higher effective vibrational quantum 
number n ~ E a /uj a w 14, than those in multiphoton excitation experiments where the energy is distributed between 
many modes. At the energies at which it is formed, SFg~ also has a stronger anharmonicity related to the proximity 
of dissociation thresholds. 

The vibronic state of SFg~ at the instant of capture, \I/att, Eq. (|12p . is embedded in the dense multimode vibrational 
spectrum. The average level spacing of the SFg" vibrational spectrum at the e~ + SFg threshold is about 10~ 10 a.u. 
(see Sec. IHIj) . In the process of IVR, the initial state \l/ a tt spreads over a large number of multiple vibrational 
excitations of SFg~. For a weak coupling (i.e., small 7), \t a tt would describe a simple single-mode vibrational Feshbach 
resonance. In the actual case of strong coupling, ^ at t plays the role of a doorway for the final multimode vibrational 
resonances. Therefore, the ultimate states of SFg" populated in electron attachment are extremely closely spaced 
metastable vibrational resonances with very large but finite lifetimes. In the next section we consider their decay via 
electron autodetachment and estimate the corresponding lifetimes. 

III. DECAY OF THE METASTABLE NEGATIVE ION 
A. Evaluation of decay widths 

The metastable negative ion species SFg"* formed by electron attachment can d ecay via the reverse, autodetachment 
process, Eq. (J4j) . Another decay channel open at low energies, see, e.g., Refs. [U HI, HH , is dissociation, Eq. ((3]). 
To analyse these possibilities, compare the dissociative attachment cross section with the measured cross section for 
the production of SFg" and calculated a a tt, Fig. [31 The relatively small SF^~ signal at £ < 150 meV (below the 
dissociation threshold) is due to electron attachment to thermally activated SFg. It depends strongly on the target 
molecule vibrational temperature and has been used as a "thermometer" [8| . At larger energies the SFJ~ cross section 
rises rapidly and at e > 0.3 eV it becomes the dominant contribution to the total electron attachment cross section. 

In our view this picture indicates that for most of its part dissociation does not occur as a decay of SFg * formed 
via the resonant doorway state \I/ a tt- If the latter were true, the calculated cross section cr at t, which is actually the 
cross section of formation of metastable SFg"*, would follow the total attachment cross section, i.e., the sum of the 
dissociative attachment cross section and that of SFg" production. We conclude that at low energies, dissociation does 
not proceed via the intermediate multimode vibrational resonances formed in s-wave electron attachment described 
by \&att- Hence, autodetachment is the main process responsible for the decay of SFg~* [48[. 

The situation when a projectile forms a very dense spectrum of resonant states with the target ("compound states") 
is well known in nuclear physics [49| . In this case it is useful to consider the cross section avera ged over an energy 
interval containing many resonances. Such average cross section is described by the optical model [401 ]. In this model 
the resonant cross section, in our case cr at t, accounts for all processes occurring via intermediate resonant states (i.e., 
SFg - *), including their contribution to elastic scattering. If the mean energy spacing between the resonances, D, is 
much larger than the resonance widths, the cross section of electron capture by the SFg target in the initial vibrational 
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FIG. 6: Ratio of the mean resonance width to the level spacing evaluated from the attachment cross section, Eq. (|20|l . for 
go = 0.034 and ki = 4.1, no = (solid curve) and no = 1 (dashed curve). 



state v can be written as 

ffatt ( e , (20) 

where T„ is the average partial elastic width at a given energy [4fj| | . It determines the average detachment rate leading 
to a free electron with energy e = k 2 /2 and neutral SFg in state v. In our approach <j att depends only on the number 
of breathing vibrational quanta in the state v, i.e., on no (see Sec. Ill A[) . Figure [5] shows the ratio T v /D obtained 
from Eq. I|20p. as a function of the electron energy for no = and 1. 

Since T v / D <C 1, the assumption of non-overlapping resonances is valid. Hence, we can evaluate D as the reciprocal 
of the level density of resonances (p = D^ 1 ), and use Eq. (|20|) to estimate the lifetimes of the metastable states, 
t = h/T. In this context Eq. (|20[) is sometimes referred to as the principle of detailed balance [l||. In addition to 
finding T^, one also needs to take into account level-to-level fluctuations of the widths, which affect the shapes of the 
SFg* decay curves (see Sec. lIIIB]) . 

Apart from the total energy, the process of electron attachment and detachment considered here conserves the total 
angular momentum J, its projection M, and parity. The electronic part -00 of the doorway state V&att, Eq. (I12p . 
is spherically symmetric, and the continuum electron is represented by the s-wave. Hence, the angular momenta of 
the SFg target and SFg~* resonances coincide. These resonances also have the same parity as the initial and final 
vibrational states of SFg. Therefore, D in Eq. (|20|) refers to the average spacing between the SFg" levels with definite 
J, M and parity. The SFg molecule is a spherical top, and its rotational states with a given J and M are 2 J + 1 
times degenerate with respect to the quantum number K [40J. For the highly-excited SFg" this degeneracy can be 
lifted by ro-vibrational interactions. The rotational energies of the anion and neutral molecule are close, and much 
smaller than the total vibrational energy of SFg~ , 

k 2 

E = E a + E n + — , (21) 

where E Vi is the target initial vibrational energy. Hence, we can write 

D- x = ±(2J + l)p(E), (22) 

where p is the total density of vibrational states of SFg" , and \ accounts for parity. 

The total vibrational energy E is large (> 1 eV), but the mean energy per vibrational mode is comparable to their 
frequencies, and we calculate the vibrational spectrum of SFg using the harmonic approximation (see Appendix iBj). 
Table U lists the mode frequencies of SFg and SFg\ In contrast to SFg, the fundamentals of SFg" are not well known. 



TABLE I: Vibrational modes and frequencies of SF6 and SF 6 . 



Mode 


Symmetry 


SF 6 ° 


Frequencies (cm 1 ) 
SF 6 " 6 


SF 6 " C 


1 


Alg 


774 


700 


626 


2 


E a 


642 


625 


447 


3 


Tin 


948 


925 


722 


4 


Tin 


616 


594 


306 


5 


Tig 


525 


500 


336 


6 




347 


325 


237 



"Experimental data from Ref. [Soil . 

'Frequencies used in Ref. [35| a nd earlier in Ref. [5lll ■ 

C MBPT(2) calculation, RefTQ. 
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FIG. 7: Cumulative level numbers of even (solid) and odd (dashed) vibrational states of SFs. The inset shows the vibrational 
spectrum density of SFe (dashed) and SF^ (dot-dashed and solid) obtained using the frequency sets from Table [I] 




In the past these frequencies were assumed to be equal to those of the neutral, or slightly softer [H, [H, HH, HE, 
For example, the first set of anion frequencies in Table [J was used in Refs. [35l l5lj|. The second set is from Ref. 
which is probably the best calculation of SFg and SF^~ by the coupled-cluster and many-body perturbation theory 
(MBPT) methods. These lower anion frequencies are supported by a recent photodetachment study of SF^ [53], and 
we regard them as more accurate. On the inset of Fig. [7] we show the total vibrational level densities in the energy 
range relevant to metastable SF^*. As expected, the density obtained using the second set of anion frequencies is 
much greater than the density given by the first set. 

The average total resonance width T is the sum of partial widths T„ for all open decay channels v allowed by the 
conservation of energy, angular momentum and parity. Equations (|20p and (|22|1 give 

T(E) = ^^f , (23) 

where the energy e v and momentum k v of the emitted electron arc determined by e„ = = fc 2 /2 + E Vi — E v . 

The sum in Eq. (|23p is over the open-channel vibrational states of neutral SFg. Summation over 2 J + 1 states with 
different K gives an additional factor 2 J + 1, which is cancelled by the same factor in Eq. (f2"2"|) . To compare with 
experiment, the width (|23|) can also be averaged over the distribution of the initial target states and incident electron 
energies. 




Electron energy (meV) 

FIG. 8: Mean resonance width T as a function of the incident electron energy e for different initial target states: ground state 
(solid curve), Ti u (dashed curve), Tig (dotted curve), and T\ u (dot-dashed line). Shown on the inset are the populations of a 
few lower vibrational states of SFg at room temperature. 



To within a factor of two related to parity selection, Eq. (f2"3"f coincides with that derived by Klots [33j using 
equilibrium considerations. It has a simple physical meaning. The decay rate in each channel is given by the 
probability per unit time for SFg~* to enter the doorway state ^attj which is equal to the classical frequency D /2n, 
times the probability P of detachment at each attempt. The latter quantity is determined by the capture cross section, 
P = CT a tt (s, v) k 2 /n, as follows from the detailed balance relations 18l |. 

To show how the number of open decay channels affects the width, we plot in Fig. [7] the cumulative number of 
SF6 vibrational excitations for both parities. The parity of the final SFg must be the same as that of its initial state. 
Besides the ground state, at room temperature only the lowest excited states, T2 U , T2 g , and T\ u , are populated with 
probabilities w u > 0.05. These are shown on the inset of Fig. [H In this figure we also show the energy dependence 
of the mean resonance width T, Eq. (f2"3")h for each of the above target states. The dependence of T on e is step- like, 
each step related to the opening of a new decay channel. The width is also bigger for higher-lying target states, as 
more channels are open. 

The widths in Fig. [8] are given in reciprocal time units, i.e., they represent the decay rates T/h. The smallest of 
the widths, for the electron capture by the ground-state SFg, is about 2 ms _1 below 65 meV. This indicates a lifetime 
of 0.5 ms, which is comparable to the values observed in traps [23|, [24|, [25|, [26|. Note however, that the widths are 
very sensitive to the SF^~ vibrational spectrum. Thus, if we used the first set of anion frequencies from Table HI the 
widths given by Eq. (|23p would be about 10 2 times larger, since they are inversely proportional to the vibrational 
level density of SFg shown in Fig. [7J The smallest width would then be about 10 /xs, i.e., in the range of values from 
time-of-flight experiments [l?], H, O HH, H3] • In fact, the earliest of such studies, Ref. [HI, which analysed the 
lifetimes theoretically in a way similar to Eq. (123|) . using SFg fundamentals and disregarding parity, concluded that 
the lifetime of 25 /is corresponded to the electron affinity of ~ 1.1 eV. 

It should also be noted that due to the strong energy dependence of p(E), the widths and lifetimes are sensitive to 
the value of the electron affinity. A change of 0.1 eV in E a changes the lifetimes by a factor of three. In addition, 
our calculation of the density neglects anharmonic effects in the anion vibrational spectrum. The amount of energy 
deposited in each mode is relatively small, but the total vibrational energy E is close to the dissociation threshold. 
In this case anharmonicity can increase p(E), and hence the lifetimes, by a factor of two or three. 

Figure [S] shows that if the electrons and SFg molecules possess broad thermal energy distributions, e.g., when the 
anions are formed in a trap, the decay of SFg~* will be characterised by a set of average widths rather than a single 
detachment width. This was indeed observed by Delmore and Appelhans 0, [2l[ , who detected several lifetimes or 
bunches of lifetimes depending on the temperature (in the 10 ps range, though) . For monoenergetic electrons with 
e < 65 meV and vibrationally cold SFg, the decay process is governed by a single average width. Such experiment 
was done by Garrec et al. [22j, who did report a single lifetime r w 19 p,s. However, the corresponding width would 
be close to the values found using the neutral or similar vibrational frequencies (35| . It is incompatible with the more 



accurate softer anion vibrational spectrum, unless a much lower electron affinity is used. 

In our view, the present calculation does explain observations of sizeable SFg~* signals at millisecond and greater 
times, and its slow nonexponential decay [24j, [25|, l26( , at least qualitatively. Equations (|20|) and ([23| yield the widths 
averaged over large numbers of closely spaced resonances. Given the large density p(E), these resonances cannot 
be resolved experimentally, even with a highly monoenergetic electron beam. As a result, one always observes an 
ensemble of such states. Some of them may have widths much smaller than the average. Such states will form a tail 
of long-lived anions, as was first suggested by Odom et al. [24]. Therefore, to determine the survival curve of SF^~*, 
one must take into account the distribution of resonance widths. 



B. Fluctuations of the widths and nonexponential decay 

The width of a particular multimode vibrational resonance is determined by the size of the doorway state \l/att 
component in its wavefunction. This component is extremely small since ty^tt spreads over a large number of such 
resonances, which can be estimated as "fp(E) ~ u> a p(E) ~ 10 10 (Sec. Ill B 2|) . Physically, this situation is similar to 
that of neutron capture by heavy nuclei, where each of the narrow compound resonances contains only a small fraction 
of the "neutron + target" state, which allows their coupling to the continuum. Due to strong mixing between the 
basis states that make up the compound states, the statistics of their components becomes Gaussian (see, e.g., Ref. 
(Hll]). This leads to the following probability density for the partial widths r„, 

e -r„/2r„ 

/(r„) = = , (24) 

where T v is the mean. Equation (|24p is known as the Porter-Thomas distribution [49|, [54| . 

If only one decay channel is open, every resonance decays exponentially as e~ Tvt . However, it is easy to see that 
fluctuations of T v result in a nonexponential time dependence of the survival probability P(t) for the ensemble. 
Assuming that at t = all levels in the ensemble have equal populations and using Eq. (|24p. one obtains, 

/>oo 

P(t)= / e- r - t f(T u )dT u = (1 + 2TA)- 1 ' 2 . (25) 
Jo 

For more than one decay channel, the survival probability is given by (see, e.g., [HBL l56l|) 

^) = n( 1 + 2T ^r 1/2 ' ( 26 ) 

V 

where the product includes all open channels. When the number of open channels, _/V c , is large, then for T v t <C 1 one 
can use (1 + 2r„i)~ 1 / 2 ~ e -1 ^* in Eq. Pa]), which gives 

P{t) ~ e- rt , (27) 

where T = r„ ~ N C T V is the total width. The exponential behaviour that holds for Tt < r/r„ ~ N c is a 
consequence of suppression of fluctuations in the total width. The case of N c ^> 1 in fact corresponds to the classical 
limit of the decay process of a compound system with a dense spectrum of states. The parameter N c determines 
whether the decay process is classical with an exponential behavior, or quantum where sizeable deviations from the 
pure exponent can be expected (see, e.g., (57j). 

Besides the fluctuations of the partial widths described by Eq. (|24p. the distribution of resonance widths in the 
negative ion ensemble depends on the conditions of its formation. Thus, SFg~* formed in a trap will have a different 
distribution of autodetachment widths and hence, different decay curves, to SF^~* formed with a high-resolution 
electron beam. In particular, the population of SF^ resonances created in a beam experiment at time t = is 
proportional to the level density and their elastic (entrance) widths. For a single decay channel v (identical to the 
entrance channel), the number of ions that have survived to time t is found as 

N(t) cx p(E) J r.e-^/Cr,)^ = (i ^ 2 ^ 3/2 . (28) 



Note that compared with Eq. (|25p , the account of the initial resonance population has resulted in an additional factor 
(1 + 2r t/ t)~ 1 . For more than one decay channel, Eq. ((2l)|) will be similarly modified, 

N(t)K^^ t l[(l + 2T„t)-^, (29) 



where Vi is the initial vibrational state of the target molecule. 

If the target molecules are characterised by a distribution of initial states with vibrational and rotational quantum 
numbers Vi and J, with probabilities w Ui j, the decay curves must be averaged over these. This gives 

« E r+ar't/^ n^ 1 + n»*m- Njl *, (30) 

where iVj is the number of channels with different rotational quantum numbers K' of the final SFg that can be 
populated for a given K of the target molecule. If the K quantum number is not conserved than Nj = 2 J + 1, while 
if K were conserved, one would have Nj = 1. In the former case the decay should be close to exponential for room 
temperature (or even much colder) SF@, since typical J and Nj are large. Indeed, the probabilities for a thermal 
ensemble of SFq molecules are w Vi j oc exp{— [E Vi + J(J + l)/2/]/T}, where T is the temperature in energy units, 
and I — 1.23 x 10 6 a.u. is the moment of inertia of SFg. Hence, from J(J + l)/27 ~ §T we have J ~ 50 at room 
temperature. 

The widths in Eq. (|3H)) depend on the initial vibrational state of the target and the incident electron energy. 
For SF^~* formed in a trap in equilibrium with thermal electrons, the resonances will be populated uniformly, with 
Boltzmann probabilities oc e~ E l T . In this case the decay curve is given by the right-hand side of Eq. (f3T)|) without 
the factor w Vi jT Vi /(l + 2T Vi t/Nj), averaged over the Boltzmann anion energy distribution. 

Figures l9l and ITOl show the results of our modelling of the decay curves of SF^~* formed under different conditions. 
In both figures the bottom plates show the apparent lifetime r estimated from the instantaneous decay rate r _1 = 
(dN/dt)/N. This quantity should be constant for a purely exponential decay, but in general, given Eqs. (|2"5|) (pfUf . 
increases with time. The rate of such increase is asymptotically linear, r oc t, since the fraction of long-lived ions 
increases with time. 

In Fig. [5] we assume that the anions are produced in collisions of monochromatic electrons with four different 
energies with the SF 6 gas at T — 300 and 10 K. As expected from the energy dependence of the widths, Fig. [SJ 
the decay becomes faster much with the increase of the electron energy. However, even at the highest energy the 
lifetime remains greater than 100 /is, which is compatible with the laser photoelectron attachment experiments 0, @]- 
Non-exponential features clearly seen in Fig. [5] at small times are due to the contribution of short-lived anions formed 
by vibrationally excited target molecules. Naturally, this effect is more pronounced at higher temperatures. The 
nonexponential decay at large times is caused by the Porter-Thomas fluctuations of the autodetachment widths. This 
effect is greater at small molecular temperatures, where the number of rotational channels is not too large. However, 
at room temperature Nj ~ 50, and the decay is close to exponential over the time interval shown. 

The curves that model the decay of SFg~* formed in a trap at two temperatures, T — 300 and 77 K, are shown in 
Fig. 1101 In addition to nonexponential features at small times, the presence of long-lived anions formed by low-energy 
electrons leads to a nonexponential decay at large times. To illustrate this effect we show for comparison in Fig. \W\ 
the decay curve for anions formed by monoenergetic electrons with the energy of |T. Also shown in Fig. [TU] is the 
decay curve calculated assuming that K is conserved (i.e., neglecting the effect of the rotational motion on the number 
of decay channels). As expected, this decay curve is strongly nonexponential, with the apparent lifetime growing as 
t ~ t. This feature of the SFq decay was observed experimentally in Refs. [H|,[2(|. It is a consequence of the fact that 
the fraction of surviving long-lived ions increases with t. On the other hand, allowing for the mixing of the rotational 
quantum number K makes for a faster (exponential) depletion of SFg~ , which makes it hard to explain how the anion 
signal can be observed at much longer times. 

There is another effect that molecular rotations can have on the SF^"* lifetimes, and that has not been taken into 
account in the present consideration. Due to a difference between the equilibrium S-F bond lengths of the neutral and 
anion, the moment of inertia of SF^~ is about 20% greater than that of the neutral. This means that the contribution 
of the rotational energy J (J + l)/27 in SFg and SF^ differs by 20% as well. At room temperatures this difference is 
close to 10 meV, which means that the amount of energy available for IVR in different rotational states is different, 
and is larger for higher J. This is another source of fluctuations in the detachment widths, which can in principle 
affect the anion decay curves. 



IV. CONCLUSIONS 



Comparison of calculated cross sections with experimental data demonstrates that low-energy electron attachment 
to SFq proceeds via capture into a virtual state with a strong coupling to the breathing vibrational mode. The model 
provides an accurate description of the attachment cross section, and is in agreement with the measured total and 
vibrational excitation cross sections. From this comparison, we have determined the two parameters of the model 
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FIG. 9: Time evolution of the survival probability of SF 6 *, N(t)/N(0) (top) and apparent lifetimes, r (bottom) calculated 
for two target temperatures T = 300 K (thick curves) and T — 10 K (thin curves), and different electron energies: 0.45 meV 
(solid), 44 meV (dashed), 79 meV (dotted), and 113 meV (dot-dashed). 

which describe the behaviour of the SFq and SF^~ potential curves near the equilibrium of the neutral. This behaviour 
is in accord with the quantum chemistry calculations of the potential curves. 

By allowing a reflected wave in the nuclear dynamics of the SFg breathing mode, we have studied the effect of 
the IVR rate on the cross sections. Comparison with experiment indicates that the IVR in SF^~ formed by electron 
attachment is very fast, its rate being comparable to the frequency of the breathing mode. 

We have evaluated the autodetachment widths of metastable SF^~* resonances, assuming statistical distribution of 
the energy over the vibrational spectrum of the molecule. The magnitude of the widths depends strongly on the set 
of SFg vibrational frequencies used, as well as on the adiabatic electron affinity of SFg. Using the recommended 
value of 1.06 eV together with the best calculated fundamentals we obtain estimates of SFg~* lifetimes in the 100 
^s to 1 ms range, shorter lifetimes corresponding to the anions formed by higher-energy electrons. These lifetimes 
are broadly in agreement with the values inferred in laser photoclcctron attachment experiments and found in ion- 
cyclotron-resonance experiments and in traps. At the same time, we cannot explain the observation of tens of /its 
lifetimes in time-of- flight experiments. Such values would be compatible with a much stiffer set of SF^ fundamentals 
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FIG. 10: Time evolution of the survival probability of SFg~*, N(t)/N(0) (top) and apparent lifetimes, r (bottom), calculated 
for anions formed at equilibrium conditions at two trap temperatures, T = 300 K (solid) and T = 77 K (dashed). Dotted curve 
corresponds to the decay of SFg"* produced by capture of electrons with fixed energy e — 38 meV, equal to the mean energy 
at room temperature. Dot-dashed curve shows r obtained assuming conservation of the angular quantum number K (i.e., for 
Nj = 1). 

(similar to the neutral), or smaller electron affinity values. 

By using the Porter-Thomas distribution we have modelled the effect of fluctuations of autodetachment widths due 
to statistical ( "chaotic" ) nature of highly-excited vibrational states of SFg~ . We have also investigated the effect of 
SFg temperature and electron energy on the SFg * decay curves. Fluctuations of the widths result in nonexponcntial 
decay of SFg~*. However, the presence of a large number of rotational channels for all but very low temperatures of 
the neutral molecule makes these effects relatively small. 

Finally, dissociative attachment into SFJ~ + F has been largely ignored in the present work. We believe that it 
can have only a very small effect on the lifetimes of metastable SFg formed at low electron energies. It appears 
that for electron energies below 0.2 eV the SF^ ions originate from (thermally activated) metastable SFg. However, 
the dissociation signal is low due to a small branching ratio in comparison with autodetachment. At higher electron 
energies the SF^ signal has a large peak. Here the dissociation cross section exceeds that of the s-wave attachment 



model. In our view this means that the main dissociation mechanism at these energies is different from that responsible 
for the long-lived SFg~ . 

When this work was completed, we became aware of two very recent studies of the electron-SFg problem. In the 
first of these [58| , the autodetachment lifetimes of metastable SFg were studied with a time-of- flight and Penning ion 
trap techniques, for very low (~ 1 meV) electron energies as a function of the gas temperature. At room temperatures 
only long-lived SFg~* with t > 1 ms were seen. At higher temperatures shorter lifetimes were observed (~ 0.4 
ms), together with a small signal of short-lived anions (< 10 /is). These values are broadly in agreement with our 
decay curve modelling. The second set of papers [59| used kinetic modelling within the framework of the statistical 
unimolecular rate theory. It analysed the attachment and dissociation data in a wide range of target and electron 
temperatures, as well as SFg and carrier gas pressures, allowing for additional effects (e.g., the rate of IVR as a 
function of the incident electron energy) by model factors and fitting parameters. One of the results of this analysis 
is an indication of a larger electron affinity of 1.20 ± 0.05 eV. This value would result in a factor of 3-5 decrease of 
our autodetachment widths, and a corresponding increase of the lifetimes. Given the sensitivity of the decay curves 
to the conditions under which SFg~* are formed, such change still leaves the lifetimes in the 1-10 ms range, for the 
conditions similar to those used in Ref. 1581 1. 
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APPENDIX A: MATCHING OF THE WAVEFUNCTION 



Here we show how the total wavefunction \l/(r, q) of the e~ + SFg system, Eq. ([H]), with the amplitudes satisfying 
Eq. (fT0|) . matches the adiabatic wavefunction ^ at t(r, q), Eq. (fT2| . of SFg~, in the range of nuclear coordinates where 
the incident electron is bound. In order to do this, we project \l/att onto the SFg breathing vibrational states Xn(q) 
with large n. This allows us to determine the asymptotic behaviour of f n corresponding to the two terms on the 
right-hand side of Eq. (|13|) and obtain an explicit relation between /„ and the amplitudes A and B. 

Consider the wavefunction \P a tt(r, q) f ar from the merging point of the neutral and anion potential curves. Here its 
nuclear part (fT3"|) can be written explicitly using the semiclassical (WKB) approximation [4Q|. Let us first consider 
the contribution of the outgoing wave term, 

= /^r^ ex P (* / P{l)dq J , (Al) 

Vv{q) V Ja J 



where p(q) — y/2M[E — Uo(q)], v(q) = p(q)/M is the corresponding velocity, E = \k 2 + E no is the energy of the 
system, and a is in the classically allowed region (its choice only affecting the phase of A). 
Let us now demonstrate that the projection of \& a tt = ^>o( p i 9)x(<?) onto Xn(q), 

9n(r)= / <Wr,«)Xn(?)dg, (A2) 



matches the electronic part of the terms in the sum over n in Eq. jS]). 

For a large n, we can use a normalised semiclassical expression for Xn(q) [sol ] . 



Xn(q) = J ^(q) C0S (/ P^^^t) (A>!l 



where p n — y/2ME n — (Muq) 2 , v n (q) = p n (q)/M, and a n = ^/2E n /Muj 2 is the classical turning point. Using Eqs. 
(|A"Tj) and in (|Al]l . we obtain: 



,(r) = ^ / &£^L exp I ; / p[q)dq 1 cos [ I Pn[q Ul q - -),/,. (A4) 



The integral above contains rapidly oscillating functions. It can be evaluated in the saddle-point approximation, and 
the main contribution is due to the incoming wave component of the cosine. Hence, the oscillatory factor in the 
integrand is of the form exp[i<p n (q)], where 



p{q)dq 



p n (q)dq. 



(A5) 



The equation for the saddle point, d<p„/dq — 0, gives p(q) = p n (q), i.e., the "transition" between the nuclear 
vibrational state of the neutral and that of the anion takes place when the momenta are equal. This gives the 
following equation for the saddle point q n : 



E-s (q n ) = E n . 



It immediately follows that 



n(q n ) = ^J-2e{q n ) = \/2(n - n )ui - k 2 = \k n \, 



(A6) 



(A7) 



which shows that g n (^) oc e l fc ™l r /r, exactly as that of the nth term in Eq. for closed channels. Note also that the 
saddle point, found using Eq. (|A3|) . as 



|fe n | — Kq \J1ion 



(A8) 



lies in the classically allowed region of the oscillator, which justifies the use of Eq. (|A3|) . 

Completing the saddle-point calculation of the integral (|A4[) , we obtain an expression for the amplitude /„ at large 
n, for the case when x(<?) is an outgoing wave: 



fn 



A 



A 



2-Kv(q n )K, 1 

2cun 
2ttki M 



exp[iip n (q n )] 



Muj* 



-1/4 



ex.p[iip n (q n )}. 



(A9) 
(A10) 



In the last equation we used an asymptotic expression for the velocity v(q n ) = \/2[E — Uo(q n )]/M obtained for large 
n using Eq. J7|. 

Thus we see that apart from a slowly varying pre-factor, the successive amplitudes differ by their phase, so that 



n+1 



In 



exp 



i-7-<Pn(q n ) 

an 



Using Eq. (|A5[) and taking into account the fact that dip n (q) / dq — at q = q n , we obtain: 



d ( ^ 
an 



a " dq 



v n (q) 



where r n is the time it takes the oscillator with energy E n to pass from q n to a n . Hence, we obtain 

q n 



£ = cos uiT n + i sin cuT n 



+ u 1 



(All) 



(A12) 



(A13) 



where 



qn 
a„ 



MW \k n \—Kg 

2n + 1 Ki 



(A14) 



We can see that expression (| A13|) coincides with ^ + from Eq. (fT5|) for large n — N . This proves that the solution of 

the recurrence relation (|10[) with the boundary condition ff^li/f^ = ^+ corresponds to the outgoing wave in the 
anion nuclear wavefunction. 

If instead of (|Alj) we consider the contribution of the incoming wave in x(?)> Eq. (113|) . the coefficients f n will 
be given by the complex conjugate of Eqs. (|A9|) . (|A10[) with A replaced by B. The corresponding phase factor 



fn+i/fn is the complex conjugate of (|A13|) . asymptotically equal to £_ = from Eq. (fT3| . Hence, the incoming 

wave contribution is obtained by solving Eq. (fTUj) with f^~h/fj^ ^ = 

The relation between the behaviour of ^ a tt(r, g) at large q and the asymptotic form of the amplitudes /„ for closed 
channels allows one to find the attachment cross section cr at t directly from /„. Thus, for the pure outgoing wave case, 
calculating the flux for x(q) from Eq. (|Alj) . and using Eq. (|A9j) . we have: 



L4f 
k 



2tt 



\fn\ 2 Kl 

ku> 



v(q n ), 



(A15) 



where the last expression is independent of n for large n. 

Finally, the reflection coefficient R is evaluated by following the semiclassical wave Eq. (|A1|) along the anion 
potential curve Uo(q) from the matching point a = q m to the right turning point b and back. Here we assume that 
due to IVR the potential also acquires a small negative imaginary part, —17/2. This gives 



R 



-1 exp 



2i / p(q)dq - 7 



dq 
v(q) 



(A16) 



where the momentum p(q) and velocity v(q) are calculated as in Eq. (|A1[) using Uo(q), and the imaginary part of 
the potential results in damping at the rate 7. If the anion potential curve for q > q m is described in the harmonic 
approximation, 



U (q) = U a 



^Mu 2 a (q-q a ) 2 , 



the integrals in Eq. (|A16|I are given by 

eh 



7 



p{q)dq = 

dq 

v{q) 



E-U a 



arcsm a 



1_ 



arcsm a 

2 



(A17) 

(A18) 
(A19) 



where a = (q m — q a )/{b — q a ), b — q a = y/ {E — U a )/(Mujl), ui a is the frequency of the SF^T breathing mode, and \U a \ 
is the adiabatic electron affinity (neglecting the zero-point energy). 

The reflection coefficient that relates the amplitudes corresponding to the outgoing and incoming waves, Rf = 
fx/fx, is obtained with the help of Eq. (|A9[) . and contains an additional phase factor: 



R 



f 



R exp 



-2i 



Pn{q)dq 



In the matching procedure the coordinate q m is chosen equal to gjv, Eq. 
the recurrence relation (fTU)) . 



(A20) 

81), where N is the truncation number of 



APPENDIX B: VIBRATIONAL SPECTRUM DENSITY 



Let us calculate p(E) as the density of states for an ensemble of s harmonic oscillators with frequencies LUi for a 
fixed temperature T, 



e s 



where 



p{E) = ; (Bl) 



S =§ -J2Ml-e- i/T ) (B2) 
»=i 



is the entropy of the ensemble, 



i—l ^ ' 



is the variance of the energy, and T is measured in the units of energy. To find the density for a given energy, the 
temperature must be chosen so that the mean energy of the ensemble (measured from the ground state) is equal to 
E: 

i—l 

The density of states calculated in this way is in excellent agreement with that obtained by directly counting the 
multimode vibrationally excited states. 
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